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Abstract. Planned space-based ultra-high-energy cosmic-ray detectors (TUS, JEM- 
EUSO and S-EUSO) are best suited for searches of global anisotropics in the 
distribution of arrival directions of cosmic-ray particles because they will be able to 
observe the full sky with a single instrument. We calculate quantitatively the strength 
of anisotropies associated with two models of the origin of the highest-energy particles: 
the extragalactic model (sources follow the distribution of galaxies in the Universe) 
and the superheavy dark-matter model (sources follow the distribution of dark matter 
in the Galactic halo). Based on the expected exposure of the experiments, we estimate 
the optimal strategy for efficient search of these effects. 



PACS number: 98.70.Sa 



Global anisotropy of cosmic rays with space-based detectors 



2 



1. Introduction 

The next step in tlie studies ultra-high-energy (UHE) cosmic rays (CRs) is related to the 
use of space-based detectors observing the fluorescent hght induced by air showers in the 
terrestrial atmosphere with large exposures. It is expected that these instruments will 
help us to shed light on the origin of the highest-energy cosmic particles which remains 
unknown up to now. 

One of the important signatures of particular UHECR models is the global 
anisotropy of arrival directions of the highest-energy events. For instance, models 
where the origin of UHECRs is attributed to the acceleration in astrophysical objects 
(so-called "bottom-up" scenarios) naturally predict that the distribution of arrival 
directions follows the distribution of these cosmic accelerators. In the most common 
scenario with extragalactic protons and/or nuclei, the patterns of the distribution of 
galaxies in the nearby Universe should be seen [1] on the UHECR skymap because of 
the hmited propagation length of these particles due to the GZK effect [2] or nuclear 
photodisintegration. Strong suppression of the cosmic-ray flux is predicted at highest 
energies in these models. On the other hand, the "top-down" models [3] (see Refs. [4] for 
reviews and references) with the distribution of sources in the Galactic halo following 
that of the dark matter (which is the case for the superheavy dark-matter (SHDM) 
particles and some of topological defects) predict [5] the Galactic center-anticenter 
asymmetry due to the non-central position of the Sun in the Galaxy (see Refs. [SJ [TJ [8] 
for extensive discussions). Models of this kind predict continuation of the cosmic-ray 
spectrum and gamma-ray dominance beyond 10^° eV. 

Currently, neither the spectrum nor anisotropy observations can deflnitely favour 
one of the two scenarios at the highest energies, E ~ 10^'' eV. Indeed, the AGASA 
experiment claims [9] the super-GZK continuation of the spectrum while the HiRes 
collaboration reports [10] the observation of the Greisen-Zatsepin-Kuzmin (GZK) cutoff 
(data of other experiments, including recent results of the Pierre Auger Observatory 
(PAO) [HI |12], are not yet conclusive: though the unsuppressed continuation of the 
spectrum is excluded by the Auger data, the cutoff in the spectrum is not clearly seen). 
The limits on the gamma-ray fraction in the primary cosmic-ray flux (currently the most 
restrictive ones arise from the AGASA and Yakutsk muon data at E > 10^° eV [13], from 
the Yakutsk muon data at E > A- 10^^ eV [H] and from the Pierre Auger Observatory 
data on the shower geometry at E > 2 ■ 10^^ eV and E > 10^^ eV [I5]) disfavour the 
SHDM scenario and even exclude it for particular values of the dark-matter parameters 
(see e.g. Ref. [16] for a more detailed discussion of some of these limits). Current 
experiments do not report any signiflcant deviations from the global isotropy at the 
highest energie^. This is however not conclusive both due to the low statistics and due 
to a limited fleld of view of any terrestrial-based installation. 

The steeply falling spectrum of cosmic rays makes it very difficult to obtain a reliable 

f After this paper was submitted, PAO reported a significant deviation from isotropy jllj. Tfie definite 
interpretation of this effect awaits further data (see e.g. Ref. [18] for discussion). 
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measure of the global anisotropy in any combination of terrestrial experiments. Indeed, 
the relative systematic difference in the energy estimation between two installations 
located in different parts of the Earth (and thus observing different parts of the sky) 
can hardly be made smaller than some ~ 15%. Such a relative error would give ~ 30% 
higher integral flux seen by one of the experiments with respect to the other at the same 
reconstructed energy. Thus, possible observations of global anisotropy can be attributed 
both to a physical effect and to unknown systematics in the energy determination. 
Moreover, a seemingly isotropic distribution of the arrival directions over the sky might 
represent a physically anisotropic one masked by the systematic effects. 

On the other hand, the planned space-based experiments, e.g. TUS [T9], JEM- 
EUSO [20] or S-EUSO |2l], will provide a unique opportunity to observe full sky 
with a single detector. While not being free from systematic uncertainties in the 
energy determination, an experiment of this kind would not introduce strong direction- 
dependent systematics and thus would be able to perform the studies of the global 
anisotropy at high confidence. 

One may expect that in future space-based experiments with their huge exposures, 
particular sources of UHECR will be determined on event-by-event basis for the case of 
the astrophysical scenario. This is not an easy task, however: limited angular resolution 
together with large numbers of events would lead to identification problems similar to 
thoose of the gamma-ray astronomjj§| but strongly enhanced due to magnetic deflections 
of the charged cosmic-ray particles. The searches for global patterns in the distribution 
of UHECR arrival directions will thus be important in any case. 

A robust method for the study of any global asymmetry in the arrival directions is 
the harmonic analysis (see e.g. Ref. [23]). It works perfectly if the predicted effect may 
be clearly seen in the first few harmonics but requires large statistics to reveal/exclude 
more complicated patterns. Here, we determine the optimal strategy for the searches 
of global anisotropies even with the low-resolution experiments (TUS) and for short 
exposure times. 

The strategy is to fix two regions of the sky (not necessary covering full An) which 
are expected to provide the strongest contrast in over/underdensity of events with 
respect to the null hypothesis of the isotropic distribution. The shape and the size 
of these regions, as well as the energy range of the events, are determined a priory in 
order to balance the strength of the expected anisotropy (increasing for smaller regions 
and higher energies) and its statistical significance. The aim of this paper is to simulate 
optimal regions and energies for TUS, JEM-EUSO and S-EUSO, suitable to test the 
scenarios of extragalactic sources and of decays of superheavy dark matter concentrated 
in the Galactic halo. To this end, we perform new and improved (with respect to 
previous works) simulations of the distribution of arrival directions expected in both 
cases. 

§ See e.g. Refs. [22| for discussions of statistical methods of analysis of photon-by-photon EGRET 
data. Currently, the gamma-ray sky at energies (~ 0.1 1) GeV contains 101 identified source, 170 
unidentified ones and a strong non-uniform diffuse background of unidentified origin. 
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The paper is organised as follows. In Sec. [2] and Sec. [21 we discuss our simulations for 
extragalactic sources and for SHDM decays, respectively. SecHlgives specific predictions 
for three coming spaceborn experiments, TUS, JEM-EUSO and S-EUSO. Sec. [5] contains 



our brief conclusions. Some technical details are collected in Appendix A 



2. The "bottom-up" scenario: sources follow the distribution of 
extragalactic luminuous matter 

A number of astrophysical sources were suggested where acceleration of cosmic rays 
up to the highest energies can take place (see Refs. [24] for reviews and summary). A 
common assumption is that the distribution of the sources follows that of luminuous 
matter, that is of galaxies. We should stress that this approach does not imply that 
all galaxies emit cosmic rays of UHE energies - this is certainly ruled out; instead, 
it assumes that the number density of the sources is, on average, proportional to the 
number density of the galaxies. This latter assumption is true for most of the suggested 
cosmic accelerators (active galactic nuclei, starburst galaxies, gamma-ray bursts etc.) 
Our approach assumes that the number of sources within the GZK sphere is large enough 
(> 10^) so that the large-scale structure of the Universe traces well their distribution 
within ~ 100 Mpc. The assumption does not work for the case of only few sources and 
very strong intergalactic magnetic fields (see Sec. 12.31) . 

Interaction of UHE hadrons with cosmic background radiation limits the 
propagation distance at high energies; hence, a limited part of the Universe may contain 
the sources of UHECRs detected at the Earth. Non-uniform distribution of matter in 
this part should reflect itself in the distribution of the arrival directions of cosmic rays [1]. 
For this study, the following scheme was used. 

• We calculated the density function of the sources n{l, b, d) as the number density 
of galaxies for a given direction (Galactic coordinates /, b) and distance d. 

• By making use of a numerical propagation code [25], we estimated the fraction 
f{d,E^i^) of surviving hadrons with E > -Emin at the distance d from the source, 
assuming either proton or iron injected primaries. 

• The density of the sources n{l,b,d) was convolved with the survival function f{d) 
to obtain the expected distribution of the UHECR arrival directions. 



2.1. The density function of the sources 

Most of the previous studies [1] used the PSCz catalog [26] to reconstruct the source 
density. However, due to a limited angular resolution of the IRAS instrument used 
to compile the catalog it might not resolve all galaxies in the regions of high density 
(clusters) [27] and thus may systematically diminish the expected anisotropy due to 
underestimation of high number densities of galaxies. For now, the most complete 
full-sky catalog of galaxies is the 2MASS XSC [28]. Redshifts are known only for 
a small fraction of the 2MASS galaxies; the photometric redshifts are proven to be 
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a good measure of the distance to the rest of them [21]. However, the precision of 
photometric redshifts is insufficient at low distances; moreover, artificial nearby sources 
may appear due to uncertainties in the determination of the Galactic extinction. We 
use another full-sky galaxy catalog, LEDA [30], to identify the density function at 
c? < 30 Mpc. A volume-limited sample of LEDA [31] for d < 50 Mpc is used, and the slice 
30 Mpc< d < 50 Mpc is used for calibration of the density of the sources (in arbitrary 
units) determined from LEDA and from 2MASS XSC photometric redshifts. The latters 
are calculated following Ref . [32] , using the Galactic extinction of Ref . [33] , and are used 
for d > 30 Mpc. The details of calculation are given in Appendix A Our sample is 
complete up to the distances of order 270 Mpc (we assume Hq = 73 km/s/Mpc). 

The actual density of the sources is much less than that of all galaxies. This 
means that the density function should be smoothed on a reasonable scale to suppress 
unphysical fiuctuations. For this study, we smooth the function n(/, b, d) at the scale 
of 10 Mpc in d and ~ 5° on the celestial sphere. The latter value is also motivated by 
the uncertainty due to defiection of charged hadrons in the cosmic magnetic fields (see 
Sec. 12.31 for discussion) and by the experimental angular resolution. 



2.2. Calculation of the fraction of surviving hadrons 

We use a detailed propagation code [25] which is based on kinematic equations written 
in the expanding Universe and accounts for numerous interaction processes, tracing the 
propagation of nuclei (iron and lighter), nucleons, gamma-rays, electrons and neutrinos. 
Parameters of the simulation have been chosen [31] in such a way as to provide the 
best fit to the HiRes data [10]: for injected protons, the assumed spectral index at 
the source is a = 2.55 and the maximal injected energy is -Emax = 1-3 ■ 10^^ eV; for 
iron these parameters are a = 2.2 and -Emax = 1-7 ■ 10^^ eV, correspondingly. In both 
cases, the best fit corresponds to the absence of the source evolution and to the zero 
distance to the nearest source. The same parameters provide reasonable fits to the data 
of the Pierre Auger surface detector [TT] if the energy calibration independent of the 
fiuorescent yield [12] is used; see Ref. [35] for other fits to the Auger spectrum. 

The functions f{d,Emi-a) are plotted in Figs. [T], [2] for different injected primaries 
and -Emin- Since our sample is complete up to ~ 270 Mpc, we see (cf. Fig. [H [2]) that 
the study makes sence for E >7 ■ 10^^ eV. 

These results are compared to a related study [36] in Figs. [3], HI The agreement 
between the two simulations is reasonably good though ours predict slightly larger 
contribution from distant sources. 



2.3. Effect of the cosmic magnetic fields 

The previous discussion did not take into account the effect of the cosmic magnetic 
fields on the propagation of charged cosmic particles. This effect is determined by 
bending of particle trajectories and is therefore twofold: firstly, the trajectories become 
longer and secondly, they do not point back to the origin. The first effect may have 
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Figure 1. The fraction of surviving hadrons with E > Emin as a function of the 
distance d from the proton-emitting source with spectral index a — 2.55 and a cutoff 
at 1.3 • 10^1 eV. Sohd hue, E^in = 4 • lO^^ eV; dotted Une, E^in = 7 ■ 10^^ eV; dashed 
Une, En,in = lO^o eV; dash-dotted hne, E^in = 2 • lO^o eV. 





Figure 2. The fraction of surviving hadrons with E > i?min as a function of the 
distance d from the iron-emitting source with spectral index a = 2.2 and a cutoff at 
1.7 ■ 10^2 eV. Solid line, E^i,, = 4 ■ lO^^ eV; dotted Hne, E^^in = 7 • lO^^ eV; dashed 
line, Sinin = lO^" eV; dash-dotted hne, i^min = 2 • 10^° eV. Note that the number of 
particles is first enhanced at dosens Megaparsecs from the source because of prompt 
disintegration of heavy nuclei into a number of light ones and only then becomes 
suppressed due to the GZK effect. 
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Figure 3. Horizons of protons: proton-injecting sources at distances larger than D 
contribute a fraction k of the observed flux. Solid lines, k = 0.3; dashed lines, k = 0.1. 
Black lines, this study; gray lines, Ref. [HQ- Injection spectral index a — 2.7. 
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important consequences for the spectrum if the change in the trajectory length becomes 
comparable to the attenuation length of the particle; the second effect may change the 
expected skymap significantly. 

For our purposes, it would be useful to separate the effects of the extragalactic 
magnetic fields from those of the magnetic field of the Milky Way. 

2.3.1. Extragalactic magnetic fields. Direct information about the strength of 
intergalactic magnetic fields is not available and the corresponding estimates vary by 
orders of magnitude. Magnetic fields of order 10^^^ G in the most part of the local 
Universe (except dense clusters of galaxies) are predicted by simulations of Ref. [37] 
while fields as strong as 10"^ G in voids are advocated in Ref. [38]. Detailed calculations 
of cosmic-ray trajectories were given e.g. in Refs. [371 EH] for both cases. A typical 
deflection angle x of particle with energy E and charge Z may be roughly estimated 
as 

7-lO^^eV B 
X^u.i ^ 10-11 G' 

where B is the typical magnetic-field strength in the voids. We see that for the strong- 
field scenario {B ~ 10~^ G), cosmic-ray trajectories are bend strongly and become 
entangled (x ~ 360°) even for protons (Z = 1) of ~ 7 ■ 10^^ eV (in this case, even 
a few sources within ~ 100 Mpc would not result in significant anisotropy of arrival 
directions). For less extreme values oi B < 10~^ G, typical deflections of protons do 
not exceed a few degrees at the super-GZK energies so that the analysis of this work 
remains applicable. For the weak-field scenario of Ref. [37] which we assume in what 
follows, deflections of nuclei as heavy as iron are of order of, or less than the experimental 
angular resolution aX E >7 ■ 10^^ eV. 

The effect of the lengthening of trajectories on the spectrum of various nuclei was 
studied in detail in Ref. |39]; it may be neglected within our precision for B < 10~^ G 
and the energies of interest (see e.g. Fig. 10 of Ref. [39]). 

2.3.2. Galactic magnetic fields. Much stronger (~ 10^® G) and somewhat better 
studied magnetic fields are present in the Milky Way (see e.g. Refs. [IQ] for reviews 
of various models with applications to UHECR deflections). However, their strength 
is partially compensated by the relatively small size they occupy (strong fields are 
localized within the Galactic disk and possibly in the central regions of the bulge). At 
E >7- 10^^ eV, the deflections are always small (~ 1°) for protons but may be important 
for heavier nuclei, at least in the directions along the Galactic plane. The actual skymap 
at low Galactic latitudes b depends strongly on the particular model of the magnetic- 
field structure; however at relatively high \b\ > 50° the deflections are always small and 
the pattern of cosmic structures at these latitudes (where e.g. the Virgo supercluster is 
located) remains unchanged. The lengthening of particle trajectories (except possibly 
for the rare ones which hit the central region of the Galaxy and may get entangled) is 
negligible as compared to the attenuation length. 



Global anisotropy of cosmic rays with space-based detectors 



9 



To summarize this section, the results of our study are not significantly affected by 
the account of intergalactic magnetic fields of 5 < 10^^ G for protons and B < 10"^'' G 
for heavier nuclei and by the account of the Galactic magnetic fields - everywhere for 
protons and at high Galactic latitudes, |6| > 50°, for heavy nuclei. The actual account 
of these magnetic fields is strongly model-dependent and is beyond the scope of this 
paper. 



3. The "top-down" scenario: sources follow the distribution of the Galactic 
dark matter 

In "top-down" models, UHECRs originate from decays of superheavy particles (the 
latters themselves may be produced in decays or collisions of the topological defects). 
The distribution of sources thus follows the distribution of the initial particles; for 
many realistic models it corresponds to the distribution of the dark matter. Decays 
of particles from the Milky- Way halo dominate the cosmic-ray flux at high energies in 
this case. At lower energies, this flux should be supplemented by a contribution of 
astrophysical sources to explain the observed spectrum. Due to a non-central position 
of the Sun in the Galaxy, the flux should be higher from the direction of the Galactic 
center than from the opposite one. 

Decays of the SHDM particles may be described in a more or less model-independent 
way because the most important physical phenomenon of relevance is hadronization 
which involves light particles and is well understood. Denote x = j^, where E 
is the energy of a decay product of the SHDM particle with mass Mx- Then for 
10"^ < X < 0.1, spectra of the decay products calculated by various methods [3111] are 
in a good agreement with each other; moreover, the shape of the spectral curve ^{x) 
does not depend on Mx |Z]. For this study, we use the spectra of decay products from 
Ref. [7jO 

The SHDM decay rate is determined by the concentration nx and lifetime tx of 
the SHDM particles, rix = nx/^x- The flux of secondary particles at the Earth is then 
determined by 

. ^. 1 dN 

Tx dh 

where 

j"'r--i^ (1) 

is the geometrical factor (r is the radius-vector from the Earth; though in principle one 
should integrate over the Universe and account for relativity effects, in most interesting 
cases the dominant contribution comes from the Galactic halo [5]). The mass Mx is 
subject to cosmological limits (see e.g. Ref. |12] and references therein); the lifetime tx 
is much less restricted. 

II We thank M. Kachelriefi for providing numerical tables of the functions caculated there. 
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Since the decays of SHDM particles may describe only the highest-energy part of 
the spectrum, the flux is assumed to be a sum of two components, one of which, F^g, 
corresponds to the "bottom-up" contribution (we use the results of Sec. [2] to model 
it) while the second one, -Fdm, is due to SHDM decays. Neglecting the contribution 
of SHDM outside the Milky- Way halo, one may rewrite the integral in Eq. ([T]) in the 
variables related to the Galaxy, namely the distance R to the Galactic center and the 
angle 6 between the directions "observer - Galactic center" and "observer - integration 
point" . This is useful because the distribution of the dark matter in the Galaxy depends 
on R. Though the actual shape of the distribution is currently a subject of debates, the 
profile 



is often used, where Rs = 45 kpc and a = 1 for a popular NFW model [13]. We take 
the size of the halo Rh = 100 kpc and the distance from the Sun to the Galactic center 
Tq = 8.5 kpc to obtain 



(note that the corresponding equation of Ref. [8] works only for < 6 < 7r/2). 

Each particular experiment has its own direction-dependent exposure A{1, b), so the 
observed flux is a convolution of Feg + -^dm with A. For instance, AGASA sees the 
Northern sky and has higher exposure towards the Galactic anticenter while Southern 
experiments (SUGAR and PAO) have higher exposures towards the Galactic center. To 
simulate the expected distribution of events, we fit the observed AGASA spectrum with 
the sum of F^g and Fqm and find the 95% CL regions for the three parameters: the 
coefficients at F^g and F^m (the latter is related to tx) and Mx- We impose additional 
constraints on the parameters. First, Mx should satisfy the bounds of Ref. [42] and 
should be at least twenty times larger than the energy corresponding to the last point in 
the spectrum being fitted (in order to always be in a safe region x < 0.1); these bounds 
are not very restrictive. However, we impose the limits on the gamma-ray fraction as 
reported in Refs. [T3l [HI [15] which exclude a significant part of the parameter space. 
Details of the fitting will be reported elsewhere; here we will present the results for two 
particular fits: the best fit for the AGASA spectrum (which is marginally consistent 
with gamma limits) and the fit with the minimal contribution from SHDM (consistent 
with gamma limits automatically) still allowed within 95% CL region of parameters. 
We note that the results for the best fit of the PAO surface-detector spectrum [11] 
(calibrated as suggested in Ref. fl2]) lay between these two. 





where 
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4. Predictions for space-based experiments 
4.I. Experiments 

The space detector TUS [19] is under construction and is planned to be launched in 
2010. The TUS lower energy threshold is estimated as 7 ■ 10^^ eV (with almost energy- 
independent sensitivity at higher energies). The UHECR arrival direction accuracy is 
different for different zenith angles 6: for 6 < 30°, the direction is estimated only roughly 
as being vertical in the error cone of 30°; for 30° ^ 6' < 60° the error approaches 10° 
and for 60° < 6* < 90° the error is less than 10°. We used the events with 30° < 6* < 90° 
in our estimates; the exposure factor per one year of operation is about 3000 km^ sr for 
these cuts. The TUS direction-dependent exposure A{1, h) is determined by the satelhte 
orbit parameters and by the zenith angle cut. It can be conveniently parametrized by 
a sum of a monopole and a quadrupole over the celestial sphere. 

The JEM-EUSO detector [20] (to be installed in 2012 on board of the International 
Space Station, ISS) is expected to have an instantaneous aperture of ~ 6 ■ 10^ km^ sr 
and a duty cicle of ~ 20%. JEM-EUSO will be sensitive to cosmic rays with energies 
E > A - 10^^ eV; the energy-dependent sensitivity is reported in Ref. pO]. The angular 
resolution of the detector is a few degrees; the ISS orbit parameters result in an almost 
uniform exposure over the sky (we therefore assume A{l,h) =const for JEM-EUSO). 

S-EUSO [21] is a proposed detector with an instantaneous aperture of ~ 2 • 
10^ km^ sr. The proposal has been submitted to the ESA "Cosmic Vision" program; the 
proposed launch date is 2017. For S-EUSO we assume the duty cycle, energy-dependent 
sensitivity and A{1, b) similar to those simulated for JEM-EUSO. 

4-2. Large-scale observables and optimisation 

To test one of the hypotheses discussed in Sec. E], [3l we divide the sky into three 
parts: part 1 corresponds to the strongest expected excess of events over the uniform 
distribution, part 2 corresponds to the strongest deficit and part 3 is the rest of the sky. 
The shape and size of the parts, as well as the energy range, are determined from a 
priory simulations in order to maximize the signal. The observable of relevance is the 
difference between numbers of events in regions 1 and 2, ni — n2 = A. Though this 
strategy is best suited for TUS with its relatively low aperture and precision, it will 
be useful for analysis of early JEM-EUSO and S-EUSO data as well, especially at the 
highest energies. 

Let Pi, i = 1,2,3, be the probability of an event to arrive from the ith region. 
Pi +P2 +P3 = 1; let be the total number of events in the sample. Then the probability 
to have ni events in the region 1 and n2 events in the region 2 is given by the multinomial 
distribution, 

N\ _ 
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E>7x10^^ eV, protons, galactic coordinates 




Figure 5. Expected density of events with E > 7 ■ 10 eV from proton-emitting 
extragalactic sources plotted on a skymap in galactic coordinates for uniform exposure. 
Darker regions correspond to larger density (linear scale). The white band corresponds 
to the zone of avoidance (|6| < 15°) cut out of the sample. Numerical values of 
the expected densities of events for this and other cases are available online from 
http : //livni . inr . ac . ru/UHECRskymaps/, 



The probability to have ni — n2 = A is 

P(A) = (i-P^-pr^'^P^'^rw^Smr^) 

X2Fi(-^,-^,l-A,(T^), 

where F and 2F1 are the Gamma function and the hypergeometric function, respectively; 
-N <A<N. 

For each particular choice of the regions 1 and 2 and of the minimal energy, we 
calculate P/v(A) as a function of N for the cases of the expected signal and of the null 
hypothesis of isotropic cosmic-ray flux (see Ref. [H] for a similar study). We expect 
that an experiment is able to confirm/exclude the hypothesis at the confidence level rj 
if 

Pn {Asign-a\ - Anuu > 0) > 1 - 7/. 

For a given rj, we determine the required number of events A^. Minimization of the time 
required to collect this number of events (estimated in each case using the spectrum 
relevant for the model we test) results in the choice of the optimal regions and energy. 

4-3. Results for the "bottom-up" scenario 

Figure \5\ presents the expected density of events for uniform exposure simulated under 
conditions described in Sec. [21 The strongest excess of events is expected along 
the Supergalactic plane, in particular from the Virgo cluster and the Perseus-Pisces 
Supercluster. Compared to the previous studies, ours suggests a stronger contrast 
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Figure 6. Expected distribution of events (arbitrary units) with E > 7 ■ 10^^ eV from 
proton-emitting extragalactic sources (thick hne) in the supergalactic latitude 6gG, 
assuming TUS exposure. The dashed line corresponds to the isotropic distribution. 

between local voids and filaments, probably due to the contrast in the source density 
function (cf. the IRAS-2MASS difference discussed in Sec. 12.11) . 

To determine the optimal choice of regions 1 and 2, we tested various cuts in galactic 
and supergalactic coordinates and their combinations. Not surprisingly, the strongest 
contrast is achieved in distribution of events in the supergalactic latitude 6sg (illustrated 
by Fig. [6]). The strongest excess of events corresponds to the band |6sg| < 15° (region 
1) while the region of the strongest deficit (region 2) is the Northern supergalactic 
hemisphere, +15° < 6sg < +90°. For any of the space-based experiments, the optimal 
minimal energy of the sample is -Emin = 7 ■ 10^^ eV, and ~ 30 events with these energies 
are needed in the full-sky sample for 95% CL evidence/exclusion of the hypothesis. 

4-4- Results for the SHDM scenario 

Not surprisingly again, the regions 1 and 2 for the SHDM search are determined as 
cones of angle ip around the Galactic center and anti-center, respectively. The optimal 
parameters are given in Table [1] for different experiments. 

5. Conclusions 

In this work, we made quantitative predictions for the global anisotropy of the UHECR 
arrival directions expected in two distinct scenarios ("top-down" and "bottom- up") of 
the origin of the highest-energy cosmic rays. 

Several refinements and improvements resulted in considerable changes in the 
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Table 1. Optimized search parameters for the SHDM scenario. The columns give: (1), 
experiment name; (2), assumed SHDM plus extragalactic fit of the AGASA spectrum 
(see Sec. |3]); (3)-(4), optimized search parameters: i?min ~ minimal energy of the 
sample, — angular radius of the cone around the Galactic Center (region 1) and 
Galactic anti-center (region 2); (5), number of events in the sample required for 95% 
CL evidence/exclusion of the hypothesis. 
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predictions as compared to previous studies. In particular, the patterns in the 
distribution of arrival directions of cosmic rays from astrophysical sources are more 
pronounced than expected before. The superheavy dark matter scenarios consistent 
with current data predict less pronounced anisotropy. 

We developed optimal observables for distinction (exclusion) of these two scenarios. 
The actual required observational time depends strongly on currently uncertain duty 
cycle, fluorescent yield and the actual spectrum. In any case, some of the scenarios 
will be tested even with the limited statistics of the first space-based UHECR detector, 
TUS, before (or at the time of) the planned launch of JEM/EUSO. However, only JEM- 
EUSO and/or S-EUSO will be able to detect the SHDM component predicted by the 
"minimal" fit of the AGASA data consistent at 95% CL with the spectrum and with 
photon limits. 

To observe, at the 95% CL, the patterns correlated with cosmic large-scale structure, 
one needs ~ 30 events in the full-sky sample with energies E > 7 ■ 10^^ eV. If 
these patterns show up at high confidence with lower statistics, this might indicate 
either significant underestimation of the particle energies or a problem in theoretical 
understanding of the origin and/or propagation of ultra- high-energy cosmic particles. 
One may try to use these patterns as a rough but independent tool for the energy 
calibration, quite important for space-based detectors. 
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Appendix A. Construction of the density function of the sources 

Appendix A.l. The 2MASS sample 

The most complete full-sky catalog of galaxies, the 2MASS XSC [28j, does not contain 
distances to objects listed there (simply because their redshifts have not been measured 
in most cases). However, one may use the "photometric reshifts" [29] assuming that 
the galaxies are standard candles in the infrared. More precisely, one assumes that all 
galaxies have the same absolute magnitudes in ii'-band, Mk = M^, (this assumption 
works good enough in average for large samples of galaxies, see discussion in Ref. j29j). 
Following Refs. [29lll5], we assume = -24.0. One has [29] 

Mk = i^corr - 5 log - k{z), 

ro 

where Kcorr is the apparent i^-magnitude corrected for the Galactic extinction, Dl{z) is 
the luminosity distance (for redshifts z < 0.1 which we are interested in, Dl{z) ^ z/ Hq] 
we use the Hubble constant Hq = 73 km-s^^-Mpc^^), tq = 10 pc and k{z) is the 
cosmic-reddening correction, k{z) ~ — 61og(l + z). The distance d to a galaxy is then 
determined as 

(i = ^ ^1 - ^1 - ^i/orolO(^— if-^ 

We use the model of Ref. [33j for correcting the observed i^-magnitude -ft'obs for the 
Galactic extinction: 

i^corr = i^obs-0.367E(5-l^), 

where E{B — V) is found for a given direction (/, b) by making use of the code 
dust_getval available from the website [53] . 

According to Ref. [29], the 2MASS XSC is a complete catalog of galaxies with 
\b\ > 5° and K < 13™. This corresponds to the photometric redshift z ^ 0.066 (that is, 
d ^ 270 Mpc). 

Appendix A. 2. The LED A sample and matching 

The photometric redshifts at low distances are not always reliable because of their 
relatively low precision (which affects the results more significantly for closeby objects). 
Moreover, the sample of nearby galaxies determined by photometric redshifts suffer from 
contamination by Galactic objects (though firm Galactic identifications are listed in a 
catalog available from the XSC website p8] and can easily be removed, one expects that a 
fraction of these objects is still present) and by galaxies for which the galactic extinction 
was determined poorly (e.g. due to molecular clouds in front of them). Therefore, one 
needs another source of data at small distances d. 

The most complete optical full-sky catalog of galaxies is the LEDA database [30] . 
It contains information about radial velocities (known mostly for nearby galaxies). The 
completeness of the catalog has been studied in Ref. pi] where volume-limited complete 
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Figure Al. The number of galaxies in our sample in slices of 10 Mpc in distance, 
as a function of distance. The dashed line represents the quadratic fit expected for a 
complete sample. 



samples were determined. We use the sample of galaxies with known distances up 
to 50 Mpc, complete for absolute 5-magnitudes Mb < —19 and \b\ > 15°. From 
this sample, we construct the source density function nLEDA(^, b, d) and compare it to 
^xsc(^; b-i d) in the region 30 Mpc< c? < 50 Mpc to determine the relative normalization 
of the two functions: 

nxsc - l.OSnLEDA- 

The normalization factor is close to unity which reflects the completeness of the both 
samples. For our study, we therefore use the function 



n(/, 6, d) 



1.03nLEDA(^ b,d), d < 30 Mpc, 

nxsc(^, ^ d), 30 Mpc < < 270 Mpc. 



Though nxsc is determined for \b\ > 5°, we have to limit ourselves to |6| > 15° because 
of the zone-of-avoidance cut in the LEDA sample. 

When the XSCz catalog [46j (see also Ref. [27]) containing distances to many 
2MASS galaxies will become available, the use of the separate database for nearby 
galaxies will be not necessary and it will be possible to include the region 5° < |6| < 15° 
in the study. 

The total number of galaxies in our sample is 212700. The completeness of our 
catalog is illustrated by Fig. lAll 
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